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Abstract 


We study a holographic p-wave superconductor model in a four dimensional 
Einstein-Maxwell-complex vector field theory with a negative cosmological con- 
stant. The complex vector field is charged under the Maxwell field. We solve 
the full coupled equations of motion of the system and find black hole solutions 
with the vector hair. The vector hairy black hole solutions are dual to a ther- 
mal state with the U(1) symmetry as well as the spatial rotational symmetry 
broken spontaneously. Depending on two parameters, the mass and charge of 
the vector field, we find a rich phase structure: zeroth order, first order and 
second order phase transitions can happen in this model. We also find “retro- 
grade condensation” in which the hairy black hole solution exists only for the 
temperatures above a critical value with the free energy much larger than the 
one of the black hole without the vector hair. We construct the phase diagram 
for this system in terms of the temperature and charge of the vector field. 
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1 Introduction 


Due to the strong/weak duality characteristic of the Anti-de Sitter /Conformal Field Theory 
correspondence (AdS/CFT) [1, 2, 3], it provides us with a powerful approach to study the 
properties of strong coupled systems by a weak coupled AdS gravity. The high temperature 
superconductivity is a potential area where the AdS/CFT correspondence is applicable. 
According to the symmetry of the spatial part of wave function of the Cooper pair, su- 
perconductors can be classified as the s-wave, p-wave, d-wave, f-wave superconductor, etc. 
From a phenomenological perspective, the onset of superconductivity is characterized by 
the condensation of a composite charged operator spontaneously breaking U(1) symmetry 
at some temperature. The holographic s-wave superconductor model was first realized in 
refs. [4, 5]. According to the AdS/CFT correspondence, in the gravity side, a Maxwell field 
and a charged scalar field are introduced to describe the U(1) symmetry and the scalar 
operator in the dual field theory side. This holographic model undergoes a phase transition 
from black hole with no hair (normal phase/conductor phase) to the case with scalar hair 
at low temperatures (superconducting phase). Holographic d-wave model was constructed 
by introducing a charged massive spin two field propagating in the bulk [6, 7, 8]. To realize 
a holographic p-wave superconductor model, one needs to introduce a charged vector field 
in the bulk as a vector order parameter. Ref. [9] presented a holographic p-wave model by 
introducing a SU(2) Yang-Mills field into the bulk, where a gauge boson generated by one 
SU(2) generator is dual to the vector order parameter. Other generalized studies based 
on this model can be found for example in refs. [10, 11, 12, 13, 14, 15]. An alternative 
holographic realization of p-wave superconductivity emerges from the condensation of a 
2-form field in the bulk [16]. 

In a recent paper [17], we have studied a holographic model by introducing a complex 
vector field p, charged under a Maxwell gauge field A, in the bulk, which is dual to a 
strongly coupled system involving a charged vector operator with a global U(1) symmetry. 
In this model there exists a non-minimal coupling between the vector field and the gauge 
field characterizing the magnetic moment of the vector field, which plays a crucial role in the 
condensate of the vector field induced by an applied magnetic field. We have studied this 
model in the probe limit at finite density. Such a setup meets the minimum requirement 


to construct a holographic p-wave superconductor model. Indeed, we have found a critical 
temperature at which the system undergoes a second order phase transition. The critical 
exponent of this transition is one half which coincides with the case in the Landau-Ginzburg 
theory. In the condensed phase, a vector operator acquires a vacuum expectation value 
breaking the U(1) symmetry as well as rotational symmetry spontaneously. Our calculation 
indicates that this condensed phase exhibits an infinite DC conductivity and a gap in the 
optical conductivity, which is very reminiscent of some characteristics known from ordinary 
superconductivity. In this sense, our model can be regarded as a holographic p-wave model. 

The probe approximation neglecting the back reaction of the matter fields is only 
justified in the limit of large q with qp, and qA, fixed. It has been shown that new 
phases can emerge (see refs. [18, 19, 20] for example) and the order of the phase transition 
can also be changed [21, 22, 23, 24, 25] once the back reaction of the matter fields on the 
geometry is taken into account. To study the complete phase diagram of our holographic 
system, we need to go beyond the probe approximation and to include the back reaction. 
While the previous paper [17] focused on the effects of the non-minimal coupling term and 
applied magnetic field on the condensate of the vector operator, in this paper we aim at 
studying the effect of the back reaction of the matter fields on the background geometry. 
We will turn off the non-minimal coupling between the vector field p, and the gauge field 
A,, since we do not discuss magnetic effect in this paper. So the model is left with two 
independent parameters, i.e., the mass m of the vector field giving the dimension of the 
dual vector operator and its charge q controlling the strength of the back reaction on the 
background geometry. We manage to construct asymptotically AdS charged black hole 
solutions with nontrivial vector hair. It turns out that depending on m? and q, our model 
exhibits a rich phase structure. 

The thermodynamic behavior of the model has a dramatic change from large m? to small 
m?. In the case with large m?, if one lowers the temperature, the normal phase will become 
unstable to developing vector hair below a critical temperature T}. The transition from the 
normal phase to the condensed phase is second order for larger q, i.e., weak strength of the 
back reaction. However, as we decrease q to a critical one, the phase transition becomes first 
order. On the other hand, for the case with small m?, no matter the value of q, there exists 
a temperature below which the condensed phase never exists. When the back reaction is 
weak, hairy solutions dominate the phase diagram below a critical temperature To through 
a second order transition, then the condensed phase terminates at a lower temperature To 
at which its free energy jumps to the one in the normal phase, indicating a zeroth order 
transition. As we strengthen the back reaction, we first encounter for a first order transition 
at temperature 7) and then a zeroth order transition at Tọ. For the sufficiently strong back 
reaction case, the condensed phase only occurs at a high temperature T > T, rather than 
at a low temperature. Furthermore, the hairy phase has higher free energy than the normal 
phase. The four critical transition temperatures T}, T2, T} and To decrease as one increases 
the strength of the back reaction. To summarize possible phases associated with different 
ranges of model parameters, we construct the phase diagram in terms of charge q and 
temperature T for a given mass. We find that the critical temperature increases with the 
charge and decreases with the mass of p,. 


This paper is organized as follows. In the next section, we introduce the holographic 
model and deduce the equations of motion of the model. In section 3, we give our ansatz 
for the hairy black hole solution corresponding to the condensed phase and specify the 
boundary conditions to be satisfied. Section 4 is devoted to calculating the free energy and 
dual stress-energy tensor. We present numerical results in section 5. For each given m?, 
we scan a wide range of q to find all possible types of phase transitions and construct the 


phase diagram. The conclusion and some discussions are included in section 6. 


2 The holographic model 


Let us introduce a complex vector field p,, with mass m and charge q, into the (3 + 1) 
dimensional Einstein- Maxwell theory with a negative cosmological constant. The complete 
action reads 


1 
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with L the AdS radius set to be unity and K? = 87G related to the gravitational constant 
in the bulk. The Maxwell field strength reads Fi, = VyAv—VvAp- Pw in (1) is defined by 
Pw = Dapo — Dopp with the covariant derivative D, = V, — iqA,. The last non-minimal 
coupling term characterizes the magnetic moment of the vector field p,, which plays an 
important role in the case with an applied magnetic field [17]. In the present study, since 
we only consider the case without external magnetic field, this term will not play any role. 

Varying the action (1), we obtain the equations of motion for matter fields 


V” Fin = iglo ol, — Pp) + tgyV" (oul, — ob oy), (2) 
D” pis = He Vy + 1970 Fup = 0, (3) 


and the equations of gravitational field 
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In the AdS/CFT correspondence, a hairy black hole with appropriate boundary con- 
ditions can be explained as a condensed phase of the dual field theory, while a black hole 
without hair is dual to an uncondensed phase (normal phase). In our case, since p, is 
charged under the U(1) gauge field, its dual operator will carry the same charge under 
this gauge symmetry and a vacuum expectation value of this operator will then trigger the 
U(1) symmetry breaking spontaneously. More precisely, we hope that this system would 
admit hairy black hole solutions at low temperatures, but no hair at high temperatures. 
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Thus, the condensate of the dual vector operator will break the U(1) symmetry as well as 
the spatial rotational symmetry since the condensate will pick out one direction as special. 
Therefore, viewing this vector field as an order parameter, the holographic model can be 
used to mimic a p-wave superconductor (superfluid) phase transition. This turns out to 
be true in the probe limit [17]: when one lowers the temperature to a certain value, the 
normal background becomes unstable and a nontrivial vector hair pẹ appears. In this pa- 
per, we continue to study this model by considering the back reaction of matter fields on 
the background geometry. 


3 Equations of motion and boundary conditions 


To construct homogeneous charged black hole solutions with vector hair, we adopt the 
following ansatz 

dr? 

—— + r*h(r)dax? + r7dy?, 

f(r) (5) 
pede’ = p(r)dz, A dz” = ġ(r)dt. 


ds? = — f (rje de? + 


We will denote the position of the horizon as r, and the conformal boundary will be at 
r — œ. Our consideration is as follows. Since we would like to study a dual theory with 
finite chemical potential or charge density accompanied by a U(1) symmetry, we turn on 


A 


A; in the bulk. We want to allow for states with a non-trivial current (Jz), for which 
we further introduce p, in the bulk. Because a non-vanishing (J,) picks out « direction 
as special, which obviously breaks the rotational symmetry in x — y plane. Therefore 
we introduce a function h(r) in the xx component of the metric in order to describe the 
anisotropy. 

The horizon r, is determined by f(r) = 0. The temperature T of the black hole is 
given by = 
fie x02 

Ar i (6) 


and the thermal entropy S is given by the Bekenstein-Hawking entropy of the black hole 


T = 


A = — rav (rn), (7) 


where A denotes the area of the horizon and Vj = f dady. 
One finds that the r component of (2) implies that the phase of pẹ, must be constant. 
Without loss of generality, we can take py to be real. Then, the independent equations of 


motion in terms of the above ansatz are deduced as follows 
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where the prime denotes the derivative with respect to r. 

The full coupled equations of motion do not admit an analytical solution with non-trivial 
Pz- Therefore, we have to solve them numerically. We will use shooting method to solve 
equations (8). In order to find the solutions for all the five functions F = {p,, ¢, f, h, x} 
one must impose suitable boundary conditions at both conformal boundary r — oo and 
the horizon r = rp. 

In order to match the asymptotical AdS boundary, the general falloff near the boundary 
r —> co behaves as 


p Pr- 5 Pes 
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where the dots stand for the higher order terms in the expansion in power of 1/r and 
yee L 1 We impose p,_ = 0 since we want the condensate to arise spontaneously. 


According ? the AdS/CFT dictionary, up to a normalization, the coefficients 4, p, Pr+ 
are regarded as chemical potential, charge density and the x component of the vacuum 
expectation of the vector operator J” in the dual field theory, respectively. 
We are interested in black hole configurations that have a regular event horizon located 
at ra. Therefore, in addition to f (ra) = 0, one must require ¢(r;,) = 0 in order for g” A, A, 
being finite at the horizon. We require the regularity conditions at the horizon r = rp, 
which means that all our functions have finite values and admit a series expansion in terms 
of (r — rp) as 
a a a ta): (10) 


By plugging the expansion (10) into (8), one can find that there are five independent 
parameters at the horizon {r}, Pern), (Tn); A(ra), X(rh)}. However, there are three useful 


'The m? has a lower bound as m? = —1/4 with Ay = A- = 1/2. In that case, there is a logarithmic 
term in the asymptotical expansion. We treat such a term as the source set to be zero to avoid the 
instability induced by this term [26]. The treatment for m? = —1/4 is very subtle. We will not discuss 
this case in this paper and instead we are going to give a detailed study in future. 


scaling symmetries in the equations of motion, which read 


eX + Nex, tot, >A, (11) 
Pz > Ap,, 2A te, Ao 7h, (12) 

and 
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where in each case is a real positive number. 

Taking advantage of above three scaling symmetries, we can first set {ra = 1, x (rn) = 
0, h(r;,) = 1} for performing numerics. After solving the coupled differential equations, we 
should use the first two symmetries again to satisfy the asymptotic conditions x(oo) = 0 
and h(co) = 1. Thus we finally have two independent parameters {p2(rp),@'(7_)} at 
hand. We shall use @’/(r;,) as the shooting parameter to match the source free condition, 
i.e., pz. = 0. After solving the set of equations, we can obtain the condensate ( J”), 
chemical potential u and charge density p by reading off the corresponding coefficients 
in (9), respectively. 

Under the third symmetry, the revelent quantities transform as 


TOAT, S3S, pod, por, pry We ua (14) 


We will use the transformation to fix the chemical potential for each solution the same, 
i.e., we work in grand canonical ensemble. 

Note that the set of equations admits an analytical solution with vanishing p,,, corre- 
sponding to the normal phase (conductor phase). This solution is just the AdS Reissner- 
Nordstrom black hole, given by 


y FO) (15) 
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with the temperature T = (3 — A and the entropy S = er? 
h 


4 Free energy and dual stress-energy tensor 


In order to determine which phase is thermodynamically favored, we should calculate the 
free energy of the system for both normal phase and condensed phase. We will work 
in grand canonical ensemble in this paper, where the chemical potential is fixed. In 
gauge/gravity duality the grand potential of the boundary thermal state is identified 
with temperature T times the on-shell bulk action in Euclidean signature. The Euclidean 
action must include the Gibbons-Hawking boundary term for a well-defined Dirichlet vari- 
ational principle and further a surface counterterm for removing divergence. Since we 


consider a stationary problem, the Euclidean action is related to the Minkowski one by a 
minus sign as 


— 2K? S Buclidean = pes —g(R a +P Lm) gE / dx V =h(2K J 4), (16) 


TCO 


where h is the determinant of the induced metric huv on the boundary, and K is the trace 
of the extrinsic curvature Kuv. ? 
Employing the equations of motion, the on-shell action reduces to 


— 2n2Se shell = 2P Vae ry TAKT — 2r — VT) (17) 


with 6 = 1/T and Vz = f dxdy. Substituting the asymptotical expansion (9) into (16), we 
obtain 
Q= T Gon shell — v2 


Euclidean D2 Ís, 


(18) 


where we have used the condition hg = y3 which can be easily found form the equations of 
motion (8). Note that for the normal phase shown in (15), one has fs = —r} — 4 any and 
h3 = x3 = 0. 

According to the AdS/CFT dictionary, the stress-energy tensor of the dual field theory 
can be calculated by [27] 


K r= 


with i, j = {t, x,y}. By using of the asymptotical expansion (9), we have 


1 
Tu = 3228s + 3h3), 
1 
K 
1 
Tyy = 3,2 (SF): 


with vanishing non-diagonal components. For the normal phase with hz; = x3 = 0, we 
find that Tys = Ty, and Q/V2 = —T,,; the former shows the isotropy in x — y plane and 
the latter gives the correct thermodynamical relation for the dual field theory to the AdS 
Reissner-Nordstrém black hole. In the condensed phase with nonzero (J,), the rotational 
symmetry is broken, thus it is expected to have Tes 4 T,,. But in both cases, the stress 
energy tensor is traceless, which is consistent with the fact that we are considering a dual 
conformal field theory at the AdS boundary. 


“In principle, we should also consider the surface counterterm for the charged vector field p,, but one 
can easily see that this term makes no contribution under the source free condition, i.e., pr- = 0. 


5 Phase transition 


In what follows we will look for condensed phases numerically. We take different m?’s into 
consideration and for each m? we scan a wide range of q which determines the strength of 
the back reaction of matter fields on the background. Our numerical results reveal that the 
system exhibits distinguished behavior depending on concrete value of m?. There exists 
a particular value of m?, for which we denote as m2. In the case with m? > m2, the 
condensed phase seems to survive even down to sufficiently low temperatures, i.e., T — 0. 
In contrast, in the case with m? < m?, the condensed phase cannot exist below a finite 
temperature. To determine the precise value for m2, we need to solve the coupled equations 
of motion (8) at very low temperatures to see whether the condensate would turn back to a 
higher temperature. The T — 0 limit is a challenge in numerical calculation. Nevertheless, 
our numerical calculation suggests that m? = 0, for which we have some to say below. We 
will consider one concrete example for both cases. In each case we find similar results for 
other values of m?. 


51 m =34 


For the case with m? > m2, we choose m? = 3/4 as a concrete example. For each value of 
q, the AdS Reissner-Nordstr6m solution always exists even down to the zero temperature 
limit. However, for sufficiently low temperature, we always find additional solutions with 
non-vanishing p, that are thermodynamically preferred. That is to say, for each value of q 
we take, there is a phase transition occurring at a certain temperature T}, where a charged 
black hole developing vector hair becomes thermodynamically favored. In the dual field 
theory side, it means that a vector operator acquires a vacuum expectation value (J,) # () 
breaking the U(1) symmetry spontaneously. Furthermore, the condensate (i) chooses a 
special direction, so the rotational symmetry in «—y plane is also destroyed. Our numerical 
calculation indicates that the order of the phase transition can be changed from second 
order to first order as one increases the strength of the back reaction. More precisely, the 
phase transition is second order for q > qe and first order for q < qe, where qe ~ 1.3575 for 
m = 3/4. 

Taking q = 1.5 > qe as a typical example, apart from the AdS Reissner-Nordström 
solution, we find another set of solutions with nonzero (Je) appearing below the critical 
temperature T}. Figure 1 presents the condensate as a function of temperature, from which 
one can see that (Jy) rises continuously from zero at T. The grand potential Q is drawn 
in the left plot of figure 2. It is clear that below the critical temperature T., the state 
with non-vanishing vector “hair” is indeed thermodynamically favored over the normal 
phase. We draw the thermal entropy S with respect to temperature in the right plot of 
figure 2. One can see that at the critical temperature Te, the entropy S is continuous but 
its derivative has a jump, indicating a second order phase transition. Our numerical results 
also suggest that the critical exponent for all q > qe is always 1/2, i.e., (Jr) ~ (1—T/T.)/. 


A qualitative change happens as we decrease q past qe. Consider the case with q = 
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Figure 1: The condensate (J,) as a function of temperature. We choose q = 1.5 and 
m? = 3/4. The condensate begins to appear at T, ~ 0.0179 and rises continuously as one 
further lowers the temperature, signaling a second order transition. 


1.2 < q. The condensate (Ia) versus temperature is presented in figure 3. Compared 
to the previous case, the condensate becomes multi-valued and we can find two new sets 


A 


of solutions with non-vanishing (J,) at temperatures lower than T ~ 1.037,, involving 
an upper-branch with large on and a down-branch with small (J,). Therefore, there 
are three states that are available to the system at some temperature, i.e., one is for 
(J;) = 0 and two for (J;) 4 0. To determine which is the physical state, we draw 
the grand potential Q in figure 4. One can find that the free energy versus temperature 
develops a characteristic “swallow tail’ which is typical in first order phase transition. The 
normal phase is thermodynamically favored at higher temperatures, but as we lower the 
temperature down to Te, the upper-branch finally dominates the system. We present the 
entropy S versus temperature in the right part of figure 4, from which one can see that S is 
also multi-valued and has a sudden jump from the normal phase to the physical condensed 
phase at T}. Clearly, the transition is first order. 

One interesting feature presented in both cases is that the hairy black hole exhibits tiny 
entropy at finite low temperatures, compared with the normal phase in figure 2 and figure 4. 
Since the values of S are obtained from the behavior of the solutions at the horizon, it is 
difficult to extract them with high accuracy at sufficiently low temperature. Nevertheless, 
our numerical results suggest that entropy remains small and smoothly decreases as the 
temperature is gradually lowered. As being a single state without any degeneracy, a super- 
conducting ground state should not have any entropy. In the gravity side, it corresponds 
to the fact that the zero temperature limit of the superconducting black holes should have 
zero horizon area, which was previously observed in refs. [28, 29]. 

The main results of this subsection are summarized by the (T, q) phase diagram shown 
in figure 5. The solid curve gives the critical temperature T, for the phase transition from 
the normal phase to the condensed phase. There is a critical value of q, denoted as qe, 
above which the phase transition is second order, while below which the transition becomes 
first order. It is also clearly that as q decreases, T, decreases gradually, which tells us that 
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Figure 2: The grand potential Q (left plot) and thermal entropy S (right plot) as a function 
of temperature. In both plots, the dashed blue curves are for the normal phase, while the 
solid curves are for the condensed phase. For T > T, one can only get the blue curve, but 
for lower temperature T < T, the condensed phase appears and has the lower free energy, 
thus is thermodynamically favored. 


the increase of the back reaction hinders the phase transition. 


5.2 m= —3/16 


Similar to the previous discussion, for each m? < mŽ we scan a wide range of q to find all 
possible types of transitions. There exist two special values of q, denoted as qa and qg with 
qa > 9g, Which divides the parameter space of q into three regions, q > qa, q8 < q < qa and 
q < qg, respectively. The thermodynamic behavior changes qualitatively in three regions. 
We may find second order transition, first order transition and zeroth order transition. 
In order to take account of the order of the phase transitions, we denote the transition 
temperature as T2, T} and Tọ, respectively. As a typical example, we will take the mass 
parameter m? = —3/16 in this subsection. In this case qa ~ 1.0175 and qg ~ 0.9537. 
Details are given as follows. 

For the case with small back reaction, i.e., q > da, we focus on the case with q = 2. The 
condensate versus temperature is exhibited in the left plot of figure 6. We immediately see 
that (J,) is multi-valued above the temperature denoted as Ty ~ 0.05972. Similar to the 
first order transition for m? = 3/4, the condensed phase has two branches, i.e., the upper- 
branch with large (J,) and a down-branch with small (J,). The free energy Q drawn in 
the right plot of figure 6 also shows a “swallow tail” shape, but it is very different from the 
one in figure 4. Comparing the free energy 2 for each solution, we find that the condensed 
solutions in the down-branch are thermodynamically favored, which only exist in a small 
range To < T < Ty. At other temperatures, including T < Tọ and T > T, it is the normal 
phase with (J,) = 0 that is thermodynamically relevant. At the temperature Tp, the free 
energy has a sudden jump from the condensed phase to the normal phase, indicating a 
zeroth order transition. The charge density p and entropy S in figure 7 show that both 


p and S are continuous but have a kink at Tə, indicating a second order transition. It is 
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Figure 3: The condensate (.J,) as a function of temperature for q = 1.2 and m? = 3/4. The 
critical temperature 7T} ~ 0.00342, is denoted by a vertical dotted line. The condensate 
becomes multi-valued at T ~ 1.037... The value of condensate has a sudden jump from 
zero to the upper part of the purple solid curve at T}, indicating a first order transition. 


interesting to note that the condensed phase terminates at a finite lower temperature Tọ. 
For the case qg < q < qa, let us consider for example the case with q = 39/40. The 
behaviors of condensate and other thermodynamical quantities are much more complicated. 
Figure 8 plots the condensate with respect to temperature, where Ue is also multi-valued 
above Tọ œ 0.03992u. But there are three sets of condensed solutions. According to 


A 


the value of condensate, we denote them as the upper-branch for large (J,), the middle- 
branch for middle (J,) and the down-branch for small (J,), respectively. We also draw 
the grand potential Q in the right plot of figure 8. We present the charge density as 
well as the thermal entropy in figure 9. The values of p and S have a sudden jump 
from the normal phase to the thermodynamically favored branch of the condensed phase 
at Tı, indicating a first order phase transition. As we lower the temperature, the phase 
with (Ja) = 0 is first thermodynamically favored, and then the middle-branch begins to 
dominate the thermodynamics through a first order transition at 7} œ 0.04102, finally the 
condensed phase ends up at the temperature To where a zeroth order transition appears. 
Our numerical results uncover that as one increases the strength of the back reaction, Tı 
decreases while Tọ also decrease but with a slowly rate, and finally T} becomes equal to 
To at qg. For the present case qg < q < qa, the value of Tı is always larger than To, so 
we have a first order transition from the normal phase to the condensed phase at higher 
temperature T} and then a zeroth order transition from the condensed phase to the normal 
phase at lower temperature Tọ. 

For the above two examples, an interesting common feature is that the thermodynam- 
ically favored hairy black hole solutions exist up to a minimal temperature Tọ, where it 
connects with an unstable condensed branch starting from higher temperature, which will 
be discussed below. In the present model we restrict the case with p, turned on, only the 
uncondensed phase can appear below To, so the free energy has a sudden jump from the 
condensed phase to the normal phase at To, indicating a zeroth order transition. Note that 
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Figure 4: The grand potential Q (left plot) and thermal entropy S (right plot) with respect 
to temperature with q = 1.2 and m? = 3/4. In both plots, dashed blue curves come from 
the normal phase, while the solid curves come from the condensed phase. Trace the physical 
curve by choosing the lowest grand potential at a fixed T. The critical temperature at which 
the condensed phase begins to be thermodynamically favored is 7T, ~ 0.003421, denoted 
by a vertical dotted line. The entropy jumps from the blue curve to the lowest branch of 
the green solid curve at T}. 


in the theory of superfluidity and superconductivity, a discontinuity of the free energy was 
discovered theoretically and an exactly solvable model for such phase transition was given 
in ref. [30]. Therefore, it is quite interesting to see whether the holographic model has 
some relation to the model in ref. [30]. 

As q decreases past qg, we see a dramatic change in the thermodynamics. We draw the 
condensate versus temperature in figure 10 for q = 19/20 and q = 9/10, from which we can 
find that hairy solutions only appear at temperatures above T;, and the general trend is that 
(J,) increases with the temperature. 3 The value of (J,) is multi-valued for the case with 
larger q. As we decrease q, this multi valuedness disappears. For completeness, we show 
the charged density p in figure 11 and the thermal entropy S in figure 12. Similar feature 
can also be found in these two figures. At a first glance, this appears to be surprising, 
since in general one expects the condensed phase to emerge at low temperatures rather 
than at high temperatures. To have a physical condensed phase, the hairy black hole 
configuration should have free energy less than the AdS Reissner-Nordstré6m black hole 
describing the normalphase. Comparing {2 between the condensed phase and normal phase, 
we can clearly see in figure 13 that the condensed phase has free energy much larger than 
the normal phase and thus is not thermodynamically favored. Therefore these hairy black 
holes represent unstable branches that do not contribute to the thermodynamics. Similar 
phenomenon was previously found in ref. [31] through a phenomenological model, known as 
“exotic hairy black holes”. Such a phenomenon also exists in some consistent truncations 
of string/M-theory in refs. [32, 33] as well as inhomogeneous black hole solutions in AdS 


3One may wonder if the curve could turn back at a higher temperature. We numerically check this 
increasing trend up to a very high temperature as we can and find the curve has a well defined asymptotic 
behavior. Therefore, we believe it will not turn around. 
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Figure 5: The phase diagram for m? = 3/4. The solid curve separates the condensed phase 
from the normal phase. The critical value qe divides the condensed phase into two parts. 
The case q > qe is associated with second order phase transition (green area), while q < qe 
corresponds to first order transition (red area). 


space dual to spatially modulated phase of a field theory at finite chemical potential [34]. 
This phenomenon of a thermodynamically subdominant condensate at higher temperature 
is known as “retrograde condensation”. 4 

We summarize the main results of this subsection by constructing the (T, q) phase 
diagram in figure 14. The upper solid curve indicates the transition from the normal phase 
to the condensed phase as we lower the temperature, which is the combination of T for 
q > qa and T; for qg < q < qa. The lower solid curve represents the minimal temperature To 
at which hairy solutions terminate. The region between the two boundary curves shrinks 
as one increases the strength of the back reaction and the two curves intersect at qg. 
The dashed curve for q < qg gives the value of Ta, above which a thermodynamically 
subdominant condensed phase appears. 

We now return to the critical mass m2. For very small m?, we can find the temperature 
To easily by directly numerical calculation. Figure 15 presents the temperature To for each 
m°. We find that the value of Tọ decreases quickly as one increases m? from its lower 
bound. But, to search for Tọ numerically becomes more and more difficult as the value 
of m? is close to the critical one. Due to the lake of numerical control at sufficiently low 
temperatures, we are not able to give the values of physical quantities, such as condensate 
(J,), charge density p and entropy S, at very low T. Nevertheless, believing Ty as a 
function of m? exhibits a well-behaved behavior and using the extrapolation, we find that 
the value of m? locates at m2 = 0 up to a numerical error °. There are also two other 


‘The terminology “retrograde condensation” was first introduced to describe the behavior of a binary 
mixture during isothermal compression above the critical temperature of the mixture [35]. A subdominant 
condensate can exist in this system in some temperature range. 

5In the case with m? = —399/160000 and q = 2, we find the turning point is Ty ~ 0.00429. For the 
case with m? = 0 and q = 2, we numerically solve the model up to the temperature T œ 0.000384 with 
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Figure 6: The condensate (J,) (left) and free energy Q (right) as a function of temperature 
for q = 2 and m? = —3/16. The dashed blue curve is from the normal phase, while the 
solid curves are from the condensed phase. To ~ 0.1487 is denoted by a vertical dashed 
line and To œ 0.05972p is denoted by a vertical dot dashed line. The condensate becomes 
multi-valued between Tọ and Ts. The lower branch of the condensed phase has the lowest 
free energy. For temperatures T < To and T > T>, the normal phase is thermodynamically 
preferred. 


hints that support our numerical result for the value of m2: (1) Comparing our model to 
the SU(2) p-wave model, the effective mass of the vector operator is equal to zero for the 
latter, we find that the reduced equations of motion look very similar with each other and 
the asymptotical expansions near the boundary are the same. Since there does not exist 
any particular temperature at which the condensate terminates and turns around to high 
temperatures in the SU(2) p-wave case, it is natural to expect that the result is also true 
in our model with m? = 0. Indeed, for the special case m? = 0, if we ignore the possible 
turning point Tọ at very low temperature, the model exhibits extremely similar behavior 
as the SU(2) p-wave model [22, 36]. Furthermore, our ongoing analysis of our model with 
the full back reaction in (4+1) dimensional black hole as well as soliton background cases 
exhibits all known phase structure presented in this paper. (2) Looking at the asymptotical 
expansion of py in (9), one can see that the leading term with coefficient p,— for m°? > 0 
is divergent as r — oo. In this case such a term must be regarded as source term and set 
to be zero. In contrast, this term for m? < 0 is regular as r — oo. Similar to the case 
in the s-wave model [4], we may have freedom to consider p,— either as source term or 
expectation value of the dual operator. The critical value of mass square is m? = 0. Thus, 
it is not surprised that the model exhibits distinguished behaviors for m? < 0 and m? > 0 
numerically. 

One may suspect that even for very large m? there exists a non zero but very small 
minimal temperature like Tọ. Indeed, we can not rule out this possibility by numerical 
approach only. Nonetheless, for large m? > m2, our numerical calculation indicates that 


A 


the physical branch of (J,) (p and S) versus temperature may behave well up to zero 


very high computational accuracy and find no turning point. 
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Figure 7: The charge density p (left) and thermal entropy S (right) versus temperature for 
q = 2 and m? = —3/16. The dashed blue curves are for the normal phase, while the solid 
curves for the condensed phase. Ty ~ 0.1487 is denoted by a vertical dashed line and 
To œ 0.0597 is denoted by a vertical dot dashed line. As Ty < T < Tə, p and S become 
multi-valued between, where only the lower branch of the charge density and the upper 
branch of the entropy are thermodynamically preferred. In other region of temperature, 
the normal phase dominates the phase diagram. Both p and S are continuous but not 
differentiable at To, characterizing a second order phase transition. 


temperature. In particular, we do not find any evidence of the branch turning back to 
a one at higher temperature as a characteristic of a zeroth order transition shown, for 
example, in figure 6. ê Of course it is helpful to clear up this issue by constructing the 
extremal limit of the hairy black hole solutions. We leave this for future investigation. 

We draw the plot of the critical temperature versus the mass of the vector field from the 
normal phase to the condensed phase in figure 15. One can see that the critical temperature 
decreases with the mass of p, at a given q. This behavior is also observed in holographic 
models involving a charged scalar field [37, 38]. 

We did not find any hint in the probe limit that the condensate as well as other 
quantities turns back to a branch extending to the high temperature region at a certain 
temperature like Tp [17]. Therefore, the appearance of To is the consequence of the back 
reaction. For sufficiently large q, the back reaction can be ignored and we should recover the 
results in the probe limit. In other words, Tọ should be zero and the critical temperatures 
T, and Ty should arrive at some finite constants as q — oo. However, one can see from 
figure 14 that the value of Tọ increases with q, and from figure 5 and figure 14 that the 
critical temperatures T, and To from the normal phase to the condensed phase also increase 
with large q linearly. This is seemingly inconsistent with the expectation in the probe limit. 
In fact our numerical results are consistent with the probe limit. The reason is as follows. 
Note that the probe limit demands that one takes the limit q — oo, while keeping qp, and 
qA,, fixed in our model. To compare our results to the ones in the probe limit, we should 


°As a typical example, we choose the model parameters m? = 3/4 and q = 43/40. We manage to 


solve the coupled equations of motion (8) up to the temperature as low as T œ 3.332 x 10~°y and find no 
evidence of the condensate turning around to a branch extending to the high temperature region. 
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Figure 8: The condensate (J,) (left) and free energy Q (right) as a function of temperature 
for q = 39/40 and m? = —3/16. The dashed blue curve is from the normal phase, while 
the solid curves are from the condensed phase. 7, œ 0.04102 is denoted as a vertical 
dashed line and Tọ ~ 0.03992; is denoted as a vertical dot dashed line. The condensate 
behaves multi-valued. The middle branch of the condensed phase has the lowest free energy 
between Tọ and T7}. For other range of temperature, the normal phase is thermodynamically 
favored. 


make the scaling transformation p, — qp, and A, — qA,. Under such a transformation, 
the chemical potential u becomes qu and the temperature T changes to T/q. Therefore, 
it is the value T/q that corresponds to the temperature in the probe limit. Indeed, in 
our numerical calculations, we checked that To/q — 0 and the critical temperatures T./q 
and T,/q approach to some constants as q — oo. Therefore our numerical results are in 
agreement with the probe limit analysis. 

Furthermore, we know from from (20) that to see whether the dual stress-energy tensor 
is isotropic or not, we have to extract the value of hs carefully in our calculations. We find 
that the numerical solutions presented in our paper have ha = 0 up to a numerical error 
(~ 1074). It means that although the model has an anisotropic structure associated with 
the p-wave order with non-vanishing (Ia) in the condensed phase, the stress-energy tensor 
is isotropic. This is the same as in the SU(2) p-wave model case ” and is consistent with 
the arguments presented in ref. [39]. 


6 Conclusion and discussions 


In this paper we studied a holographic p-wave superconductor model in a four dimensional 
Einstein-Maxwell-complex vector field theory with a negative cosmological constant. The 
complex vector field p, is charged under the Maxwell field. Taking the back reaction of 
matter fields into consideration, we managed to construct hairy black hole solutions which 
satisfy all asymptotic conditions. We found the model presents a rich phase structure 


"The anisotropy of the stress-energy tensor in the SU(2) model is controlled by the constant f? appearing 


in equation (21) of ref. [22]. It has been confirmed that the value of f? should be vanishing up to a 
reasonable numerical error [39]. 
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Figure 9: The charge density p (left) and thermal entropy S (right) versus temperature for 
q = 39/40 and m? = —3/16. The dashed blue curves correspond to the normal phase, while 
the solid curves to the condensed phase. 7, ~ 0.04102y is denoted by a vertical dotted 
line and To ~ 0.03992: is denoted by a vertical dot dashed line. When To < T < T, p and 
S are both multi-valued. In both plots, it is the middle branch that is thermodynamically 
preferred. In other region of temperature, the normal phase is thermodynamically favored. 
Both p and S are not continuous at T}, but rather jump from the blue dashed line to the 
middle branch of the solid line, signaling a first order transition. 


controlled by the mass m and charge q of the vector field p,. We investigated possible 
phase transitions in detail. It turns out that there exist zeroth order, first order and second 
order phase transitions in this model. Hairy black holes were also found in the unusual 
higher temperature range T > Ta, which always have free energy higher than the normal 
phase. The phase diagrams in terms of the temperature and charge were constructed. 

Our numerical calculation suggests the existence of a critical m? denoted as m? = 0. 
When m? > m?, we have a second order phase transition from the normal phase to the 
condensed phase for the weak back reaction case. This transition becomes a first order one 
as we increase the strength of the back reaction. The transition temperature T, decreases 
as we decrease the value of q, which means that the increase of the back reaction makes 
the transition more difficult. When m? < m?, the thermodynamic behavior of the system 
changes a lot. Starting from the high temperature region, one can find the following 
transitions: For q > qa, the system undergoes a second order phase transition from the 
normal phase to the condensed phase at Tə and as the temperature decreases to To, there 
is a zeroth order transition back to the normal phase; For qg < q < qa, the system first 
undergoes a first order phase transition from the normal phase to the condensed phase 
at Tı, then at the lower temperature To, it comes back to the normal phase by a zeroth 
order transition; For q < qg, we can only get hairy black hole solutions that are always 
subdominant in the free energy referred to as “retrograde condensation”. Here the concrete 
values of gq and gg depend on the mass squared m? of the vector field p,. 

It was argued in ref. [40] that the holographic free energy can be thought of as a sort 
of generalized version of Landau-Ginzburg free energy. In Landau-Ginzburg theory, it is 
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Figure 10: The condensate (J,) as a function of temperature for q = 19/20 (left) and 
q = 9/10 (right). The condensate only emerges above the temperature Tp ~ 0.03766 
for q = 19/20 and T, œ 0.03174 for q = 9/10. The condensate in the left plot behaves 
multi-valued. 
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Figure 11: The charge density p as a function of temperature for q = 19/20 (left) and 


q = 9/10 (right). The dashed blue curves come from the normal phase, while the solid 
curves from the condensed phase. 


usually assumed that the quadratic term depends on the temperature linearly while the 
fourth order term is not strongly temperature dependent. In our present study, we found 
the behaviors deviating from the mean field theory. In the context of Landau-Ginzburg 
theory, such deviation would be a sign of an unusual temperature dependence of the higher 
order terms. It should be stressed that our model is dual to a strongly coupled system. 
A priori the dual system does not obey the usual assumption for the free energy. It is 
in principle possible that there is some new temperature scale at which the coefficients of 
higher order terms change their signs, giving rise to non-standard phase transitions in the 
framework of Landau-Ginzburg theory. 

Our study can be straightforwardly generalized to the higher dimensional case and 
other gravitational backgrounds, such as the AdS soliton backgrounds which can mimic 
the superconductor /insulator phase transition [41]. In a recent paper [42], we studied the 
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Figure 12: The thermal entropy S as a function of temperature for q = 19/20 (left) and 
q = 9/10 (right). The dashed blue curves come from the normal phase, while the solid 
curves from the condensed phase. 
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Figure 13: The free energy Q as a function of temperature for q = 19/20 (left) and q = 9/10 
(right). The dashed blue curves come from the normal phase, while the solid red curves 
from the condensed phase. The free energy of the condensed phase in the left plot forms 
a typical swallow tail. In both plots, the condensed phase has the free energy larger than 
the normal phase, and thus is not thermodynamically preferred. 


effect of an applied magnetic field effect on the AdS soliton background, and found that the 
magnetic field can induce the AdS soliton instability due to the non-minimal coupling of 
the vector field and the background magnetic field. By comparing our complex vector field 
model to the SU(2) p-wave model with a constant non-Abelian magnetic field, we found 
that the SU(2) p-wave model can be recovered by the restriction m? = 0 and y = 1 in our 
model with the ansatz in ref. [42]. It suggests that in some sense, the charged vector model 
is a generalization of the SU(2) p-wave model to the case with a general mass squared 
m? and gyromagnetic ratio y for the vector field. Due to the adjustable parameter m, we 
can see in this paper that our model shows a much richer phase structure than the SU(2) 
p-wave model, thus can be used to describe more phenomena in dual strongly coupled 
systems. 

In the present paper, we limited ourselves to a simple case with p, non-vanishing only. 
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Figure 14: The (T, q) phase diagram for m? = —3/16. The green region is related to the 
condensed phase through a second order transition from the normal phase, while the red 
region is associated with the case by a first order transition. For other areas, the normal 
phase is thermodynamically favored. The transition temperature T, for the retrograde 
condensation is denoted by the dashed line in the left down corner. 


In principle, in order to understand the full phase structure of the model at fixed chemical 
potential, one should search for the dominant thermodynamic configuration not only in 
this given sector but in a more general setup, especially turning on the temporal component 
of the charged vector p;. This would of course be much more involved, since one should 
search for the hairy black hole configuration with the least free energy among all possible 
configurations. We will leave this issue for further study. 

As a phenomenological approach, we consider the model as a p-wave superconducting 
(superfluid) one. Indeed, this toy model could be applicable in a wide variety of condensed 
matter systems and beyond. Indeed, as we have discussed in ref. [17], it may also be 
revelent for holographically mimicking the phenomenon that the QCD vacuum undergoes 
a phase transition to an exotic phase with charged p-meson condensed in a sufficiently 
strong magnetic field [43, 44]. 

As we have mentioned above, according to the symmetry of the macroscopic wave 
function or condensate of Cooper pairs in the real superconducting materials, the super- 
conductor can be classified by s-wave, p-wave, d-wave and so on. The holographic s-wave 
model has well studied (especially with back reaction) in the literature. Adopting the 
present p-wave model, it is quite interesting to study holographic models with multiple su- 
perconducting order parameters, including the competition or coexistence between s-wave 
order and p-wave order or between two p-wave orders. 8 We will leave all these issues for 


It might be difficult to combine the s-wave model [4] and the SU(2) p-wave model [9] in one theory, 
since the order parameter in the SU(2) model is charged under a U(1) subgroup of the Yang-Mills field and 
this U(1) subgroup can not play the role of the gauge group for the s-wave model in a natural manner. The 
current p-wave model is charged under a U(1) gauge group which also can be naturally taken as the gauge 
group of the s-wave order. The competition and coexistence between two s-wave orders were first studied 
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Figure 15: The transition temperature from the normal phase to the condensed phase 
(upper curve) and the turning temperature Tọ (lower curve) with respect to m?. Here 


q = 2. The leftmost two points correspond to the mass square slightly above m? = —1/4. 


further study. 

Finally we like to mention that with suitable parameters m and q, our model shows 
the normal/superconducting/normal phase transition (see figure 14) as one lowers the 
temperature continuously. Such a phase transition is called reentrant phase transition 
in the literature [47]. The reentrant phase transition usually happens in the binary and 
multicomponent liquid mixtures. But it is interesting to note that such a phase transition 
also appears in some superconducting materials, for example, granular BaPbo.zsBio.2503 
compound [48] and cuprate superconductors [49]. Thus it would be of some interest to see 
whether our model is relevant to these superconducting phase transitions. 
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